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Using the general framework of quantum field kinetics we consider new principles to compute 
initial distribution of quarks and gluons after the first hard interaction of heavy ions. We start by 
rewriting the integral equations of QCD in the form which is generalizations of the familiar QCD 
evolution equations. These equations describe both space-time- and (s, Q^)-evolution before the 
collision, and allow one to use the ep DIS data without reference to parton phenomenology. New 
technique generate perturbation theory that avoid double count of the processes, does not contain 
an artificial factorization scale, and does not require low-momentum cut-offs since infrared behavior 

- - - is controlled by the DIS data. 
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O ' 1. Introduction. 

. Recently, there have been many calculations to find the initial distribution of quarks and gluons in heavy 

' ion collisions, which may then lead to the creation of a plasma. Contrary to other types of matter under extreme 
fSJ conditions, the situation in heavy ion physics is unique, since the very existence of quark-gluon plasma (QGP) is 
inseparable from the process of creation of the matter it consists of. Any subsequent elements of the behavior, which 
\ can be described in classical language, such as thermalization or a hydrodynamic regime for the QGP, like for any 
■ other many-particle system, is possible only after intermixing of the initial phase-space distribution. However the 
' process of creation of quarks and gluons is essentially quantum by its origin. In the present paper only this process 
, is considered. 

' Any strict formulation of a quantum-mechanical problem requires an exact definition of two main elements: the 
T-H initial state of the system, and the observables in the expected final state. This reminder would hardly be necessary 
if the main difficulties were not associated with these points: On the one hand, it is unknown how stable nuclei of the 
[ initial state are build up from quarks and gluons; on the other, there is no clear understanding of what the final state 
t may look like. So even the formulation of the problem unavoidably contains some uncertainty which requires special 
rS i' care in defining these elements. Let us begin with a discussion about what kind of problems we may anticipate. 
JL 1. The initial state: Ideally, as in, for example, an atomic collision, we would define the initial state via its wave 

function. The wave functions of QCD-nuclei are unknown. A natural alternative is to use the density matrix to 
(~| , describe each nucleus. 

I!" ' Here we may assume nuclei to be well shaped objects. The uncertainty of their boundaries does not exceed the 
typical Yukawa interaction range. In the laboratory frame both nuclei are Lorentz contracted up to a longitudinal size 
Rail ^ 0.1/m. The tail of the Yukawa potential is contracted in the same proportion. The world lines of the nuclei 
5-H ' are two opposite generatrices of the light cone that has its vertex at the interaction point. No interaction between 
nuclei is possible before they overlap geometrically. For this reason the total density matrix of two nuclei is a direct 
product of the two individual density matrices. The spaces of states where they act do not overlap either. 

Since no exact information about the initial state of the nuclei is available, it seems reasonable to rely upon the 
following two considerations: First, detailed information is inessential as it basically relates to the interactions which 
maintain every nucleus as a QCD bound state. The energy of the collision is incomparably higher. It is thus enough 
to require that density matrix yields the given total momentum and baryonic charge as averages of the corresponding 
field operators. Second, the residual dynamical information must reveal itself in the same way as in other inelastic 
processes at extreme energies, like deep inelastic electron-proton or muon-nucleus scattering. This statement may 
appear trivial, because structure functions of DIS are always used for account of this information. However, one 
should keep in mind that their definition - which does not refer to the parton model - is valid only for the DIS process 
itself. In order to apply structure functions to other interactions, using the parton model is considered unavoidable. 

The first priority of this study is to avoid any intermediate phenomenology. We insist that any information taken 
from parallel experiment is valuable only as long as both processes can be described by the same theory and with the 
same initial data. 

2. The final state: The final state is assumed to be some distribution of free quarks and gluons in the perturbative 
vacuum. This vacuum is considered to be the true ground state, and is free of QCD-condensates. It is a product of 
the nuclear collision and is postulated to exist in a sufficiently large volume. The spectrum of possible states forms 
the continuum, which is unoccupied at the beginning of the collision. 
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It is supposed also that the information which is most important for understanding the future evolution is concen- 
trated in the single-particle distributions of quarks and gluons. These distributions must be calculated from their 
quantum-mechanical definitions, keeping in mind how they are to be measured in a hypothetical experiment. Two- 
and more-particle distributions should be defined as independent elements corresponding to other measurements. 

In this work we rely heavily on a previous paper |^ , where integral equations of QCD were derived without assuming 
that averaging is performed over a stationary state. To some extent these equations resemble the diagram technique 
of Keldysh which was designed for non-equilibrium processes. They are of the same matrix form, but are not 
derived with a view to obtaining quasi-classical kinetic equations. Fields and their correlators remain the main objects 
of these equations and no phase-space distributions is introduced. For this reason the approach was named "quantum 
field kinetics" (QFK). 

These new equations are the result of an initial rcsummation of the perturbation series for probabilities of inclusive 
processes (or any other observables) . On the one hand, these new equations create an approach which proves to be 
very effective for studying various inclusive processes. On the other, it allows one to trace the temporal evolution of a 
colliding system, beginning from preparation in the past right up to the moment of measurement. This latter feature 
makes the new approach extremely attractive for our goals: The evolution of any physical system is completely defined 
by the initial data, and the act of measurement only selects one or more of all the possible quantum trajectories. 
Therefore, one can expect that results of similar types of measurement will be similar. In this paper, we want to 
consider two similar processes, viz., deep inelastic electron-proton scattering and "deep inelastic pp- or AA-coUisions," 
in parallel. 

To begin with, we must define the observables of these processes in the same way. Definitions of inclusive amplitudes 
and inclusive probabilities to find one quark or one gluon in the final state are given in Section 2. The DIS cross- 
section, viewed as an inclusive process where nothing is measured except the momentum of the scattered electron, is 
defined in the same terms at the beginning of Section 3. The rest of this section considers an instructive example of 
the lowest order calculations using a model density matrix. This calculation allows one to introduce all spinor and 
vector functions without extra complicated notations, and to clarify the complete calculation. 

The dynamical equations in their full tensor-spinor form are derived in Section 4. It immediately becomes clear 
that these equations have a ladder structure. It appears that the well known ordering of ladder cells by Feynman x 
and virtuality is a direct consequence of the retarded temporal ordering inherent to these equations in their coordinate 
form. Smaller x and bigger virtualities correspond to the later times. Thus, as a by-product, we obtain an answer 
on a very old question about such correspondence [0-^. 

The new equations appear to be richer than usual QCD evolution equations for DIS structure functions, as derived 
from the renormalization group approach. The new equations interconnect two invariant functions of the vector field, 
and two of the quark field. The new evolution equations do not depend upon the type of last interaction, however 
they may be projected onto any definite process. Specific properties of the electromagnetic interaction select only 
one of the spinor functions into the definition of the structure function F2 of unpolarized DIS. However, both quark 
field functions remain in the evolution equations, along with the two functions of the vector field. The relative scale 
of additional terms and their possible role is examined in Appendix 1. It is also shown that among other extensions, 
the new equations naturally include the BFKL | [To[ | equation and the effects of quark and gluon shadowing at small 
X. Corresponding shadowing terms appear to be parametrically larger than in previous derivations ||ll] , [l2|] . 

The meaning of the objects that obey the new evolution equation comes to light in a discussion of renormalization. 
The renormalization group approach can not be used here, if only because most of the terms in these equations 
correspond to observables (imaginary part of self-energies) are finite and may not be renormalized. Instead, we 
renormalize the second subgroup of the evolution equations for the real part of the self-energies using a conventional 
BHPZ scheme. Ultraviolet divergencies are compensated for by counter-terms from the original Lagrangian. The 
running coupling appears precisely from the requirement of renormalizability. For the moment, this part of the study 
is at the level of a basic idea. 

Using some structural, rather than quantitative, assumptions, and after projecting onto specific observables of the 
e-p DIS, the new equations can be reduced to the system of CLAP equations |l3| ~ |l5| . The objects which enter the 
new equations are similar to self-energies, and we shall call them sources. It appears that the QCD evolution proceeds 
in a way such as to create a source of a field which interacts with the detector in a certain way. The evolution causes 
the dynamical assembly of a special wave packet which represents a bare quark or gluon. This process takes place in 
real time and ends at the moment of interaction. 

We expect DIS to provide the dynamical information about this process. This information is valuable only as long 
as no measurements were done before the last interaction. All unobserved information (and providing it was not 
observed) is included in the definition of the sources with their full dependence upon x and Q^. 

Here, we do not adhere to the picture of "wee partons," and do not share the opinion that the QCD evolution 
equations describe how one valence quark develops a cloud of virtual quarks and gluons around it at small distances. 
We believe that they give (perturbatively) the "evolution" of the detector response, provided the trigger includes the 
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requirement of a bare on-shell quark in the final state (the resonant condition which results in xbj = xp). That 
they definitely do not correspond to the evolution of a single quark, is seen, for example, from possibility of including 
fusion. 

The difference between the structure functions and the sources can be explained using an analogy: in condensed 
matter or scattering theory we introduce two different quantities: the density of states, p{E), and the number of 

states below energy E, n{E) = p{E)dE. The structure functions then correspond to n{E), while the sources 
correspond to p{E). In an experiment we measure n{E) (which is proportional to the allowed volume in the phase 
space), rather than p{E). From this point of view it is not surprising that we eventually express observables, roughly 
speaking, via dG{x,Q^)/dQ^ . Different measurements study these quantities, integrating them with their specific 
"upper limits." The boundary conditions for x- and Q^-evolution are imposed by the measurement, rather than the 
initial conditions which are controlled by momentum conservation, sum rules, etc. These must be set in the past 
without any relation to the (a:;-(5)-evolution. 

Cross-sections of inclusive single-quark and single-gluon production in the lowest nonvanishing order are calculated 
in Section 5. All of them contain scale-dependent and -independent terms. The latter are much bigger than the 
former, and as a result the cross-sections are expected to exhibit only weak scale dependence. The lowest order 
cross-sections are strongly peaked at low rapidities and low transverse momenta. 

The next order of the perturbative expansion generated by the new evolution equations is examined in Section 6. 
The new expansion does not lead to any diagrams which duplicate those already included in the definition of the 
sources (or structure functions). Any such diagrams would carry severe collinear singularities. We carefully examine 
the infrared finiteness of the diagrams that do occur. It is shown that they are infrared safe and that no artificial 
cut-offs are necessary to find the total cross-section. Higher order perturbation terms are not expected to present any 
difficulty, as their final state infrared behavior will be shielded by the final distributions themselves. 

We are led to the approach advocated in this paper in an unavoidable manner. Ideally, one would start with a 
complete relativistic quantum description of the static proton and its interaction with the detector. Unfortunately, 
the many attempts to describe the bound states of QCD (see Refs. |jl^ and for the reviews) have not yet met 
with real success. To calculate the production of particles in hadronic collisions, one is limited to reasoning along the 
following lines: first, an OPE-analysis of the DIS data, which gives the structure functions of DIS; next, a partonic 
interpretation of the structure functions; and, lastly, using the factorization technique. In the end, one still faces 
severe theoretical problems caused by the soft processes, the arbitrariness of the factorization scale, etc. Here, we try 
to avoid these problems "experimentally," by maximizing the use of dynamic information hidden in the DIS data. 



2. Single-particle distributions of the partons. 

Let two nuclei A and B, with momenta Pa and Pg, move towards each other at almost the speed of light, and let 
the center-of-mass system coincide with the laboratory frame. We assume that the center-of-mass energy is very large, 
s >> Af^, so that the laboratory frame is the infinite momentum frame for both nuclei. We intend to compute one- 
particle distributions of quarks and gluons created after the first interaction of these nuclei. The inclusive amplitudes 
leading to creation of one quark or one gluon from the initial state \in) are as follows, 

{X\d{p,a,i)S\in), and {X\c{k, X, a)S\in) , (2.1) 

where d^{p,a,i) is a creation operator for an on-mass-shell quark with momentum p, spin a and color i. Similarly, 
the operator c^(k. A, a) creates an on-mass-shell gluon with momentum fc, polarization A and color a. Summing the 
squared moduli of these amplitudes over a complete set of uncontrolled states \X), and averaging over the initial 
ensemble, we find inclusive spectra of quarks and gluons 

^=Y,Spp^nSU\p,<J,i)d{p,<J,i)S , (2.2) 



^ -;^5pp,„5tct(k,A,a)c(k,A,a)5 . (2.3) 

A, a 

The initial state of the cohiding system consists of two Lorentz-contr acted nuclei which are causally independent, 
and thus the total density matrix is a direct product of two independent density matrices, 

^ PA<E) PB <E)\Q cont cont 

\. (2.4) 
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Matrix elements of pi„ are obtained by sandwiching it between all state vectors \in) which enter definition (2.1) of 
inclusive amplitude. The density matrices pA and pB contain only bound states of quarks and gluons in the presence 
of vacuum condensates. The latter are assumed to be destroyed in the course of the initial hard collision and replaced 
by the perturbative QCD vacuum. Initially, all states in the continuum are unoccupied. It means that pm contains a 
projector \0cont){0cont\ on the vacuum state in the continuum. So we may commute the quark Fock operators with S 
and S*^ and only commutators survive in the final result: 

In this expression ip^^J - {x) is the Dirac wave function of a quark. Summing over spin and color we get 

^-Ey^^-A(^^Tr[^.EJ^,(.,,)], (2.6) 
where the full 2x2 matrix of the quark self-energy is given by |^ 

S^b(x, y) = i(-l)^+^.g2 (-1)''^'' / d^dvt^rGARix, O^'disit '?)D|a,a^(^, ^) • (2-7) 



This formula implies that both quark and gluon correlators, Gar and D^.'J^ are averaged with the density matrix 
Pin given by Eq. (2.4). In the first approximation we may replace the exact qqg vertex by the bare one. Then Eq. ( |2.6| ) 
takes the following simple form. 



n dN„ g2 r ^aj^ 



dpd^x 2(27r)3 J {2tt 



{Tr[{^ + mrrGif{p-k)t'rT)['^!';{-k)] + U) ^ (B)]} (2.8) 



where the additional superscript (A) or (B) denotes that the correlation function is averaged over the initial state of 
nucleus ^ or _B, respectively. 

A similar procedure yields the following expression for the inclusive gluon production: 

- E/ <'*-"''»^4".l-*'i-.n;->,„)i (2.9) 

where the primary definition of the gluon polarization tensor Hqi is given by 
Its general formula was derived in |^: 

n>x'si^,y) = *(-i)-4+v E(-i)''^''[-/^^^'?^''G^«(^'«s,B(e,^;y)GsA(?7,x) + 

R,S=0 

+ J d^dvV^T/i^,U)T>''/n"'i^^a'^fs%tc'f'iC,^ , (2.11) 

and we postpone its further expansion because of the complexity of the emerging polarization structure. However, the 
main idea remains the same as for quark production: in first approximation we get a product of two quark, Gp'f^G^^^ 
or two gluon, Dq'^-'D^^'' correlation functions. Each of them is averaged with the density matrix of only one of the 
two nuclei. This is in line with the independence of the initial states of the colliding nuclei. 

To continue with the straightforward calculations, we should now specify an explicit form of the nuclear density 
matrix. To do this, we should solve the confinement problem - which is not our intention here. Moreover, we 
have already argued that most of the detailed information is unnecessary. Now we shall motivate this point by a 
"quantum kinematic" analysis of the extreme case when no dynamical information is required to obtain a qualitative 
understanding of the phenomenon. 
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Let the nuclei collide at an energy of about 100 TeV per nucleon, when the nuclear longitudinal size is only 10~^ fm. 
This size is much less than any known in nuclear interactions, and one may consider the domain where the nuclei 
overlap as a plane surface (or the point in the {t, z)-plane). All the subsequent dynamics takes place within the future 
light cone of this point, — > 0, t > 0. As the translation invariance in t- and z~directions is manifestly broken 
by the initial conditions, one should look for the appropriate quantum numbers (other than and p^) to describe 
the final states of the particles. The symmetry that does survive is Lorentz invariance, and the boost defined by the 
operator 

is an good quantum number. The corresponding wave functions of the free scalar particles obey the Klein-Gordon 
equation, 

I d dy^ 1 d^i, , 

(r — ) --— + m^ip = , (2.13) 



and are of the following form: 



r dr dr 



^•^^P^ (^) = • (2-14) 

Here, we use coordinates ( t — TCOsh(f>, z — rsinh0, r±), and denote = + p'±- The eigenbmctions £,^^pj^{x) 
are normalized on the space-like hypersurface r = const within the future light cone of the collision point: 

f *d 

{C,p^:^i^',p'J = / Td(t>(frC,p^i^)i-g:^Ci^',p'Jx) = v')5{p±_ -p^. (2.15) 

The boost of the particle is a legitimate, but seldomly used quantum number. It is more common to use the 
momentum. Let us try to find the functions which would behave as free plane waves, at least asymptotically. The 
wave packets with the required behavior are 

= e-^^'UpA^) = ^e-™--™^'^(^-^)e^^^--^-. (2.16) 

These wave packets represent plane waves confined to within the future light cone of the collision point, and can be 
rewritten as follows: 

Se,p^(a;) = ^ I -rp"t+rp'z^tp^r^ ^ ^ ^Qgj^^^ ^ ^2.17) 

At large m±T, the phase of the wave function Se.p^ is stationary in a very narrow interval around (p = 6 (outside 
of this interval, the function reveals oscillations with exponentially increasing frequency): the wave function describes 
a particle with rapidity 9. However, for m^r ^ 1, the phase of the wave function is almost constant along the surface 
T = const . The smaller r, the more uniformly the particle is spread along the light cone. Up to distortions caused 
by the finite size of the interaction domain, any high-energy collision will produce a distribution which is uniform in 
rapidity in the vicinity of the light cone. The picture looks as if the incoming nuclei carry this distribution ah initio. 
The latter is not surprising as the same arguments can be applied to the states of particles before a strongly localized 
interaction. The distribution dN ~ const x dd corresponds to dN const x dxp jx p in terms of the Feynman variable 
xp. Thus, we arrive at a result which is typical for the Williams- Weiszacker approach. The full consideration for the 
QCD-nucleus has been recently given by McLerran and Venugopalan Q. 

We assume that deviation from this ideal distribution can be studied perturbatively with increasing accuracy the 
more inelastic the collision is. We expect that the desired information about this deviation can be obtained from 
the deep inelastic electron-proton scattering data. We argue that these data do not imply exact knowledge of the 
density matrix of a nucleus. In order to incorporate information obtained from DIS we begin by computing the DIS 
cross-section in terms of the quantum fields kinetic (QFK) js). 



3. Deep inelastic scattering on the electron 
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The goal of this section is to find those elements of a theory which are common to two essentially different problems, 
viz.^ nucleus-nucleus (or proton-proton) collisions and deep inelastic electron-proton scattering. Our primary demand 
is that these elements should appear as a by-product of the two independent lines of calculations, initiated separately 
from first principles. 

We divide this section into two parts. For the sake of completeness we begin with a brief definition of the DIS 
cross-section in terms of the QFK-approach and define the null-plane variables which will be used for all following 
calculations. 

Before turning to a detailed derivation of the self-consistent equations in the next Section, we give an instructive 
example of the lowest order calculations. These have no direct physical value, but they allow one to overcome technical 
problems and avoid premature discussion of the highly nontrivial approximations. 



3.1. Basic definitions for DIS. 

As was emphasized in the Introduction, it is important to have similar definitions of observables for all processes 
which will participate in the future information exchange. We may rewrite Eq.( p.l[ ) for the inclusive amplitude of 
DIS as 

{X\aik')Sa\k)\tn) , (3.1) 

where k and k' are the laboratory frame momenta of the electron before and after the scattering. If g = fc — A;' is the 
space-like momentum transfer, then the DIS cross-section is given by 

da _ la L^,{k,k')W'^-^{q) 
°dk' " (4^)2 (fcP) (q2)2 ' ^■^■^) 

where W'^'^{q) is a standard Bjorken notation for the correlator of two electromagnetic currents, 

^"''('?) = ^^^H<o(9)] • (3.3) 



We accept without any discussion its standard tensor decomposition 

2xBj 2xBjM^ ' 



vWr vWri 



where i/ = qP, Q'^ = —q^ > 0, xsj = /2v and 

. ^,^^_g,.+ ^ 9 +1 ^ -q^f-L- . (3.5) 
1/ 

Hereafter we will perform all computations using the infinite momentum frame fixed by the null- plane vector n'^, 

= (l,Ot,-l), n2 = 0. (3.6) 

It defines the -components of the Lorentz vectors, 

na = = = ci^ ^ = = — 

In the infinite momentum frame, the 4-vector of the proton's momentum has components 

= (P+/2,Ot,P+/2), p-=p0-p3^0. (3.7) 

The momentum transfer has the components 

g^ = (i./P+,qt,-i'/P+), g+ = 0, q-=2iy/P+. (3.8) 

Instead of the invariant W2, we will use the mass- independent structure function F2{xBj, Q^) = i'W2/M^, which is 
calculated via the equation 

C2 = W^Uf^n, - ^ (3.9) 
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The longitudinal structure function Fl {x bj , Q2) = Wl should be calculated in accordance with 

3Fl = 2xBjCi + 2F2; ci = W^g^,. (3.10) 

3.2. An instructive example: the low order calculation. 

In order to proceed with the calculations we must specify the density matrix. We shall begin the physical motivation 
of our choice by reminding the reader that a widely used approach based on the Wilson's operator product expansion 
(OPE) docs not utilize any information about the proton's internal structure. Only the total momentum and the 
discrete quantum numbers are controlled by sum rules. Indeed, the dynamical equations of QCD contribute only to 
the singular coefficient functions while the regular operator functions, the averages over the proton's state, remain 
unknown. We can only decide whether or not to include the high twist operators in the expansion. The twist-one 
operators of OPE correspond to one-particle matrix elements of the proton's density matrix. Including the twist-two 
operators into OPE would correspond to irreducible two-particle correlations in the density matrix used here. 

The system of integral equations which we expect to derive eventually does not require any explicit form of the 
density matrix either. Nevertheless it is useful to keep in mind some representation which may serve as a simple 
reference point. For example, we may chose an artificial exponential form which reproduces the total momentum flux 
of the proton and allows to derive the integral equations of the Schwinger-Dyson type |^ . 

The twist-one operator functions of OPE, by their structure, are binary products of quark and gluon fields and to 
some extent resemble occupation numbers which enter the on-mass-shell Greenians. For example, in the statistical 
ensemble we usually have 

Gio(p) = ~2m{^ + m)6{p^ - m^Mpo){l - n(+\p)) - e{-po)n^-\p)]. (3.11) 

We define matrix elements of our one-particle (twist-one) density matrix by a certain set of the Greenians. We 
assign the superscript "#" to all states in the continmmr of free on-mass-shell fields: 

G*:'{p) = -2m5,,{^ + m)e{±p")S{p^ - m"), (3.12) 
01 

Df-^4^-(p) ^ -27ndabd^-{p)e{±po)S{p^). (3.13) 
01 

These states are initially empty and the vacuum correlators G^q2(p) and -D^Qj^(p) represent only on-mass-shell 
particles in the final states. 

The superscript will label "bounded" states of "valence" quarks and gluons in the initial proton: 



.(27r)3l 1 

■ Oij — 

Vlab 3 2 



G*oi(p) = 2m^-S,,- ^Sip')6{pt)9ip+)V{p+), (3.14) 



G*io(p) = 0. 



(3.15) 



Equation (3.14) describes the p henomenological distribution of the "valence" quarks as a function of their light-cone 
momenta p^ , while Eq. (|1|) means that ther e are no "valence" anti-quarks within the proton. (This assumes 
that we start at very low scale. Otherwise Eq. ( 3.15 ) will have the same form as Eq. ( 3.14 ).) The factors 1/2 and 
1/3 correspond to averaging of the distribution over the quark spin and color, respectively. The factor {27T)^/Viab 
corresponds to normalization: we consider a flux with one proton in a volume Viab per unit time. 
In the same way we define the "initial" distribution of "valence" gluons by 



D 



ip) 



^2ml^^Sabl:d^''{p)0{±Po)S{p^)S{pt)0{p+)g{p+), 

o Vlab ^ 



with 1/8 and 1/2 standing for the color and polarization average and where d'"' is a projector 



(3.16) 



(3.17) 



which is a sum over the physical gluon polarizations in the null-plane gauge. 
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n^B" = 0; = 0. 



(3.18) 



The reason for introducing these distributions is to give definite values to the quantum numbers (the charges and 
the momenta). Their densities are given by 



p+ 



j + {p) = Tn+G*;'{p) f dp+V{p+) = ^ f dxV{x) 

Vlab Jo Vlab Jq 



for the quark's hght-cone charge flux, and by 



T++ = -i 

I (2.) 



/ i%Tn^P^G*r{p) = 7^ r dp+p+v{p+) = ^ / 

J [^T^r ^lab Jo Vlab J a 



dxxV{x) 



(3.19) 



(3.20) 



for the (H — h)-component of the quark energy-momentum tensor. In these equations we have introduced the Feynman 
variable, x = xp ~ p'^ /P^ . The momentum flux density from the gluon component is given by 



7^++ 
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■ I ^ip^)'9,'^D;r{p) r dp+p+g{p+) = ^ f dxxg{x) (3.21) 

J (27r)^ Vlab Jo Vlab Jo 



where P^V{p'^) = V{x) and P'^Gij)'^) = g{x). The initial quark and gluon distributions are normalized in such a 
way that in aggregate they carry the proton's total quantum numbers. 



Neglecting any corrections to the electromagnetic vertex we may rewrite Eq. (3.2) in the following way, 

d^p 



^ 47r 



(2^)' 



-TrrGw{p + q)YGoiip) 



(3.22) 



The off-diagonal quark Greenians (field correlators) in this equation obey integral equations g which express them 
in terms of exact retarded and advanced propagators and sources Sio;oi (the "current" correlators): 



Gio — GretG ii^GioG ii^Gadv ~ Gret'^iaGadv ■ 

01 01 01 

The retarded and advanced Green's functions obey more familiar equations. 



(3.23) 



(3.24) 



adv 



G adv Gadv^adv^a' 



dv 



which allows a symbolic solution 



G^„, — G„g( Sret, ^adv ~ ndn, ^adv 



(3.25) 



(3.26) 



In the first approximation we may replace the exact retarded and advanced quark Green's functions, Gret and Gadv, 
by the bare ones that carry the same leading light-cone singularity, 



Gret{p) 



(pO + iQf 



ip) 



Then in virtue of ( 3.15 ), Eqs. ( 3.23 ) take the following form 

01 — ^01 ' "-^Ol ^ '-^ ret ■^01'^ adv , 
Gio = Cf^ — Gret^loGadv 

In the lowest order we neglect sources and leave only correlators of initial fields: 
W^-^{q)^ef-yi^ I ^Tn^G*{p + qh'^G*oM 



(3.27) 



(3.28) 
(3.29) 



(3.30) 



Substituting Eqs. (3.13) and (3.14) and rewriting the delta-function of the on-mass-shcU final-state quark as 
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we get the incredibly simple result of the "naive" parton model: 

/■I 

F2^\xBj) = / dxS{x — XBj)xV{x). 
Jo 



(3.31) 



(3.32) 



In the next approximation we sh ould i nclud e the quark fields coming from the quark and anti-quark sources, Sqi 
and Sioj i.e., the last terms in Eqs. ( ^3.29 ) and (3.29). The general expression for the quark self-energy matrix is given 
by Eq. (2/7). Still restricting ourselves to bare vertices and bare tree Greenians we get: 



J:**ip)^^g'CF J ^{rr7^GSi(fc)7.D#'''^(fc-p)+rr7^<(fc)7.D-^(fc-p), 



(3.33) 



(3.34) 



The superscript means that one of the contributing states belongs to the set of final states, while the other 

originates from the initial proton. 

The gluon correlators obey the following equations: 



# 



01 10 ^ Dretlloi^loDacit), 



(3.35) 



D,et — Dret + Dret^ret^ret, 



(3.36) 



"^adv — -^adv H~ -^adv^adv^ai 



■dv ■ 



(3.37) 



In the first approximation, all contributions of the gluon sources, IIio^oi, should be dropped along with radiative 
corrections to the retarded and advanced gluon propagators. The tensor structure of the self-energy may be only of 
the following form 



S(p) ^^a2{p)+ M'iip)- 
In all the equations it enters in a single combination, 

^S(p) i) ^^cri{p)+ jh.p'^p+asip); ai{p) = p^(T2{p) + 2{p+ f(73{p). 
After a long, but routine calculation wc find that 



2a** (p) = - 



VlabP+ 



Jp+ k+ p+ l~ z 



2z^ -2z + l 



(3.38) 
(3.39) 

P+g{k+)}, (3.40) 



ia**ip) = 



VlabP+ 



/ -T:^S[ip+ -k+)p' -p^,]—{CFP+v{k+) + {i-z)p+g{k+)}, (3.41) 



where z = p'^ /k^, and the invariants for the anti-quark source Sio do not contain terms with valence distributions 
V{k+). 

It is now straightforward to find the first correction F2^\xBj, Q^) to the DIS structure function: 



F!^^\xBj)^e) I dx5{x-XBj)xA^^\f{x). 



It is presented in the same form as the zero-order term, ( ^.32 ), with Vf{x) = (lf{x, Qg) replaced by 



(27r)3 



[p 



,212 



(3.42) 



(3.43) 
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To be consistent with the resonant condition of the measurement (3.31), we must require that is large enough 
(formally, — >_oo), and that the beha vior o f the integrand at high guarantees the convergence of the integral. If 
we substitute (3.41) into the r.h.s. of (3.43) and perform a residual integration over p~ using the delta-function: 



21 



- P^) '""' '^T^, (3.44) 



then we easily recover the first approximation of the Altarelli-Parisi equation for the non-singlet quark structure 
functions of the deep inelastic electron-pr oton scattering. 

In the next order we must iterate Eq. (3.35| ), including the influence of the source IIq^ on the gluon field. Cutting 
the accuracy of calculations in Eq. (2.11) to bare vertices, and neglecting the sources in the internal Greenians, we 
get the result 

Kt^'^v) = ^al{\ ^Tr[^^Gl,{k)^^GUk-p)+rG*{k+pYrGl,{k)] - 

(04K'^c/'(P> k - p, -k)-Dlf^AkKSf' i~P,P- k, k)D*^}fik - p)]} , (3.45) 

where a sum over the quark flavor / is assumed in the first term. 

The polarization tensor 11^" only appears between retarded and advanced propagators: [Dret{p)^{p)Dadv{p)]'^'^ ■ 
The latter contain projectors d^'^{p) which are orthogonal to the 4- vector rt^. So of the general tensor, only two terms 
survive: 

n^-(p) = g'^-'wiip) + ^W2{p). (3.46) 

Others, like p^^n'^ + n^p'^ or n^^n'^ , will cancel out. Introducing one more projector, 

d^^p) = -d^''ip)d;ip) ^ + ^ -P'^^ (3-47) 

which is orthogonal to both vectors n"^ and p**, we find that the invariants wi and W2 can be found from two 
convolutions, 

- d^,(p)n^''(p) = 2wi(p); n^n^W^^ip) = ^^W2{p), (3.48) 

independently of other invariants accompanying the missing tensor structures. The new projector, which includes 
only two transversal gluon modes, naturally appears in the tensor with a gluon source : 

[dip)nip)dip)r ^~dMw,{p)+^^W2{p)n''n-. (3.49) 
Now it is easy to find the first approximation for the invariants wi (p) and 'W2 (p) ■ 

x{Ci.i^-^^^-^P+V(fc+) + 27V,[z(l ~z) + ^ + i— i]p+g(fc+)}, (3.50) 



iw**{p)^--^p^ / ^5[[p+~k+)p--pi] — {ACp P+V{k+) + m,z{l--)P+g{k+)}. (3.51) 

VlabP^ Jp+ k+ p+ z 2 

In accordance with the previous convention, and as a reminder of the approximations involved, the invariants carry 
superscripts . These indicate that the invariants are contributed to by one proton's "bound" state and one on- 
mass-shell final state of a quark or gluon. We hope that the reader is not confused by the absence of other indices 
like quark color, or indices indicating the type of ordering in the invariants Wi and ct^. They can easily be recovered 



10 



when it is needed. For example, Wi and cr^ are parts of the self-energy which are summed over the color. So if 11 and 
S appear as internal elements in any formula, we must restore the color factors in the following way: 



Now we may reconstruct a missing element, viz., the first correction to the gluon structure function. 



dpt / dp p 



[p 



212 



(3.52) 



which is similar to the correction ( |3.44| ) to the quark structure function. In a sequence of approximations it should 
be added to the "valen ce" g luon distri bution G{x). 

If we substitute Eq. (3.50) into ( 3.52| ), we immediately obtain the lowest order approximation of the second Altarelli- 
Parisi equation for the gluon structure function of the deep inelastic electron-proton scattering. 

Concluding this section, let us pay special attention to the infrared poles at z = 1, which originate from the pinch- 
poles of the gluon Greenians in the null-plane gauge. To cure this problem we will proceed following Altarelli and 
Parisi ||l3]. We will first shield the IR singularity by introduction the "plus-distributions," 



^ 1 fiz)dz ^ I' /(z)^/( l) 

(1 - Z)+ Jo 



and then modify the end-point behavior in such a way that the first integrals could not be changed by the radiative 
corrections. The first integrals are the total flux of the flavor /, 



-iV, 



lab 



d^p 

(27r)4 



Trj+[Gl,{p)-G{„{p)] 



f 



and the total flux of the light-cone momentum, 

T++ = -^^a, J ^{^Tr7+p+[G^i(p) + G{„(p)] + (p+)VD(;r(f)} 



(3.53) 



(3.54) 



As follows from Eqs. (3.43) and (3.52), the first radiative corrections to the flavor and momentum fluxes are 



P+Vlab 

(27r)3 



dp^ 



dpi I dp p^ 



i(T('°^{p) i(Ji'^"{p) 



[p2 



[p2 



(3.55) 



and 



AT++ = 



(27r)3 



dp+l ^ dpi ldp-{p+f 

Qq 



E 



.iai^^{p) ia('^'^{p). iwi^{p) 



[p2 



[p2 



(3.56) 



respectively. The superscript "#*" is omitted because these equations remain valid beyond the first order approx- 
imation. The resulting conditions, Ajj" = and AT++ = 0, for the splitting kernels in the leading logarithmic 
approximation are obvious, 



Pqq{z)dz = 0, 

[zPgg{z) -f 2nfzPqg{z)]dz = 
1 



(3.57) 



[zPgq{z) + zPqq{z)]dz = Q. 

However, beyond the LLA they may change. One now readily finds an explicit form of the splitting kernels. 
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(3.58) 



zil-z) + 



where the factor /?□ = 11 — 27iy/3 coincides with the first coefficient of the Gell-Mann-Low function. 

4. Evolution of the sources. 

Now we are ready to derive integral equations that govern the field correlators and their sources. Actually, they 
have already been given above. An examination of the calculations in the previous Section shows that we did not 
sum any series. We were consequently performing a series expansion of the previously derived self-consistent solution 
of the integral Schwi nger- D yson equati ons . These are Eqs. ( p^ ) and ( 2.11| ), which define the self-energies via the 
Greenians, and Eqs. (B.29), (3.29) and (3.35), which define the Greenians (field correlators) via the self-energies. We 
shall rewrite the last equations in the form 



Gio,01 — ClfgQi — GretSio^olGa, 



dv 



and 



D 



10,01 



D 



# 



(4.1) 
(4.2) 



omitting the ^-labelled terms as they do not contribute to the differential form of evolution equations. In contrast to 
the integral evolution equations, the differential equations do not require any information about the initial data. At 
this point we do not make any approximations. 

4-1- Dynamical equations in the leading logarithmic approximation. 

In order to obtain the equations of the leading logarithmic approximation, wh ich su m up the perturbation series 

and ( ^.11 ) as the bare ones. We 



with the leading logarithms, we must consider the vertex operators in Eqs. (2 



must also confine one of the off-diagonal field correlators to the out-states in the continuum : 
^oiip) ^ ^9^CFJ-^TT{J^Gret{molik)Gadvikh.D*''^^ (4.3) 



(2^ 



TT[rGretik)J:oi{k)Gadv{kh''G*{k -p)+ -i^'G*^{k+p)YGret{k)^m{k)Gadv{k)] + 

d^k 



^ - P' ^) [Dret{k)n,,{k)D,,,{kX:, Vj;,(-P,p - fc, -k)Dt,^}f{k - p)} (4.4) 

By inspection, these equations reveal an astonishing result - the equations which govern the dynamics of the sources 
Sgi and Ilgi of the field correlators Gqi and Dgi have a ladder structure. This result appeared though we did not 
try to set any momentum or angular ordering of the emission processes. 

This result deserves spe cial discus sion . First, let us trace through the previous calculations once more. We started 
from the basic definitions (3.1) and (3.2) of the observable cross-section of the DIS. Then we expressed a specific off- 
diagonal polarization tensor of a proton via the off-diagonal Greenians Gio.oi and Dio,oi of the "kinetic" technique. 
Afterwards we used the Schwinger-Dyson equations for these Greenians to express them via their sources Sio.oi and 
IIio.oi- This resulted in a closed system of ladder-type equations for these sources. 

Now we must answer two questions: (i) what was the physical input; and (ii) what follows from the application of 
new method to the well-known problem of DIS, for which the solution is already known. 

As in any quantum-mechanic problem, the physical picture is specified by a density matrix of the initial state. Our 
density matrix carries the same information as the operator functions of the OPE, that is, no detailed information. 
Only the global quantum numbers are under its control. 

A more careful analysis shows that our result is not exactly the same as that of the OPE-based approach. Though 
we get (up to a quantitatively inessential shift of the singularities) the leading logarithm approximation resulting 
in the Gribov-Lipatov-Altarelli-Parisi equations, which sums the leading logarithms of the perturbation series, an 
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important qualitative difference appears. The Feynman tree propagators were replaced by the retarded ones. This 
change reveals the causal structure of the whole process. The last interaction, which puts a single quark onto its mass 
shell in the perturbative vacuum, is also the latest in time. In other words, the last interaction results in a collapse 
of the initial wave function. 

So, as a by-product, we have answered an old question |^-|^ about the correspondence between the evolutionary 
scale or x and the temporal scale. The causal space-time structure is contained, a priori, in evolution equations 
like the GLAP equations. It is not necessary to impose it a posteriori. The x-ordering is a further consequence of the 
0-functions that allow only for the emission into the initially unpopulated continuum. The Q^-ordering with respect 
to transverse momentum is not a necessary condition, and we shall discuss it shortly. 

What are the practical consequences of this picture? Firstly, in order to find the structure functions of the deep 
inelastic scattering of a proton off an electron, one should know the intensity of the quark field source Sqi before the 
field interaction with the electron. That is why it does not matter before what kind of interaction. The intensities of 
the sources S and 11 turn out to be universal functions. Unlike the structure functions of DIS they do not depend 
from the particular choice of the measurement procedure. We can continue this reasoning in application to pp- and 
AA-collisions. This leads to the conclusion that if we wish to use e-p DIS data for the description of p-p collision 
dynamics, we should rely on the sources S and 11, rather than the structure functions q{x,Q^) and G{x,Q^). Their 
evolution is only a specific projection of the more complicated dynamic of the sources. 

The second consequence is that after the structure functions qix, Q^) and G{x, Q^) of the e-p DIS are found (simply 
by fitting data, for example), we do not need their phenomenological interpretation as a parton density in order to 
apply them to other types of collision. 

It is well known that instead of calculating the observable e-p cross-section, the OPE method computes the imaginary 
part of the truncated Feynman amplitude of the auxiliary Compton process. The next step is a renormalization group 
analysis of this definite S-matrix amplitude. Unfortunately, the power of this method is restricted to one single 
problem. Any extension of this method requires a parton language. 

Indeed, though the cross-section of the Drell-Yan process is defined by the same (except for kinematic region) 
polarization operator as in DIS, it can not be calculated via OPE. The difference is that now the operator functions 
should be averaged over a state with two protons. Feynman propagators which contribute to an auxiliary S-matrix 
amplitude do not disappear outside the light cone. For the massless partons, and especially in the infinite momentum 
frame, it leads to the effective interaction before the collision. The factorization theorem pq | may be a remedy, but it 
requires parton language. At the same time it is clear that protons colliding at high energies are causally independent 
until the moment of collision. 



We shall now show that the above ladder equations (O) and (iA) are equivalent to the well known QCD evolution 
equations. Indeed, let us rewrite Eqs. (4.3) and (4.4) in terms of their tensor components. Beyond the first order 
calculations of Section 2.3. we must also take into account radiative corrections to the retarded and advanced 
Green's functions. The spinor (tensor) structure of Sret.adu (p^ret ad y), a s it is given by Eqs. ( |3.38 ) and ( 3.46D , remains 
unchanged. The solution of the Schwinger-Dyson equations ( 3.24 ), ( 3.25 ) and ( 3.36 ), ( |3.37 ) for retarded and advanced 
Green's functions is easily cast in the form: 



G 



ret^adv 



R.A 



ip) 



'^ret,adv^P') ^ 



/ip-^cr^'^{p) 



9 
P'^ 



R,A 
1 



R,A 



ip))' 



9 -R.A/ 
p- -W^ 



\p) {P+?p^-w^^^{p) 

Introducing the following shorthand notations for the denominators of the propagators of different modes, 



(4.5) 



W'^^\p)=p'- 



S^-^{P)^P' 



R.A 



ip), 
ip), 



s. 



R,A, 



2 {P)=P -W. 

R,A, X -I 
2 (P) = 1 - '^2 



2 

R.A 



R,A( N 
ip), 



(4.6) 



we easily obtain 



Gret (p)Soi {p)Gadv (p) 



MpV^p+)]<{p) , 



s?ip)sHp) 



^01 



ip) 



2p- 



s^siipY 



(4.7) 



[Dret{k)T\oi{k)Dadv{k)]'''' = 



di"'{p)w'(^{p) {p^l{p+Y)wf(j))ni'n'' 



W^ip)Wf{p) 



m\p)m{p) 



(4.8) 



The complete evolution equations are long, and are given in Appendix 1 along with analysis of further approxima- 
tions. Here we write only the leading logarithmic terms which eventually result in the GLAP equations: 
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-r-r 
p+ 



^2 rP^ ^1.+ r 



(2- 



fc+aOi(fc) 



p+ k+wf{k) 



k+'W^{k)Wt{k) 



(4.9) 



(27r)3 



j ^ —J d^ktdk-d[{p+ - k+){p- - k-) - (p* - 



k+a'l'{k) 
k+'S^{k)Sf{k) 



p+. k+w1\k) 



k+'W(Hk)Wfik) 



(4.10) 



Substituting Eqs. (4.9) and (4.10) into Eqns. ( ^?9| ) and ( p. 22 ), we easily find that the weU known struct ure fu nction s 
of the electron-proton deep inelastic scattering are defined by the equations which are similar to ( 3.43 ) and ( 3. 52 ): 



VlabP+ 

(27r)3 



Ql 



dp] 



dp 



ip'^a^ (p) 

WipWiP)' 



(4.11) 



G(x,( 



G{x,Ql) 



VlabP+ 



dPt I dp 



ip^Wi{p) 

w^ip)wHp)' 



(4.12) 



The first of these equations is exact, it does not depend on any further approximations. The second one is approximate. 
It holds only in the LLA. 

We observe that the structure functions of deep inelastic e-p scattering in the LLA require only one of the two 
invariants from each of the sources 11 and E. This nice feature fails beyond LLA, and does not hold for processes with 
different polarization properties of the interaction vertex. 

Unfortunately, our knowledge of the confinement dynamics is not sufficient to proceed without additional qualitative 
physical input. We must integrate over the null-plane momentum p~ and this is a highly nontrivial procedure, unless 
we postulate an ad hoc requirement that the answer should contain only logarithms. When we were doing similar 
calculations in the lowest order, we designed the density matrix of "valence" quarks and gluons in such a way that 
k~ = Q and kt = 0. In doing this, we kept in mind both the physics of the infinite momentum frame and the proton's 
confinement before the collision. 

The point is that propagation of quark and gluon fields before they reach a collision vertex is not free. Assuming 
the opposite, we would immediately violate the causality principle, or be in contradiction with previous calculations. 
We shall see shortly that the requirement of selfconsistency between emission- absorbtion and propagation provides a 
natural condition for renormalization. Suppression of the spectral patterns with fc~ 7^ is a part of this condition. 
Physically, it means that the field correlators like n(a:: — y) do not depend on the difference x'^ — y'^ of the coordinates 
in the direction of the light-cone propagation. We will conjecture that the exact retarded and advanced propagators 
indeed prevent the wave packet of a proton (or a nucleus) from a premature decay, and shall cast the requirement in 
two forms. 

The first, weak form, of this condition can be cast in the form of the inequality: \k^\ « \p^\. Then we may 
integrate over using the previous formula (3.44) as a physical approximation. The second, strong form replaces 
the inequality by the exact condition p^ — Q, which can be incorporated into the prescription: 



ip' 



ip) 



S^ip)S^ip) VlabP- 



-5{p-) 



dqf{x,Pt] 
dpi 



(4.13) 



ip'^w\ {p) 



Wf(p)Wi^(p) VlabP+ 



5{p 



_^d G{x,pf, 
dp} 



(4.14) 



These equations are a recipe on how to use DIS data in the leading logarithmic approximation. They require an 
additional comment about the pj-integration, which is the next step in obtaining the DIS structure functions. As we 
have already argued, an unambiguous definition of the structure functions is possible only in the limit of , 



00 



which eventually leads to the resonant condition ( 3.42| ) of the measurement : xp — xsj ■ 

When the r.h.s. of the Eqs. (|4.9D and Eq. (4.10) are integrated over p] with a large upper limit Q^, the kernel of the 



resulting equation depends on fcj only in the combination kf + {k^ /p'^)Q^ . The fc?-behavior of the structure function 
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in the new integral should guarantee its convergence. Therefore, one may neglect fc^ from the very beginning by 
assuming that only the domain kt <^ pt contributes in the the initial equations. This condition is known as "ordering 
by angles." Unlike the ordering by Feynman a;, it is not a fundamental r equi rement . Cle arly, this approximation is 
not valid at very low x. Nevertheless, even for low x the initial equations (4.£) and ( 4.10 ) remain unchanged. In the 
Appendix 1, we show that at low x they may be reduced to the BFKL equation. 



An explicit expression for the longitudinal structure function of the DIS follows from Eq. ( p.lO[) 



3Q^Fl{x,Q^) 



(27r)3 



dp^t 



dp 



CT2(p) 



p^Mp) 

S^{p)Sfip) S^{p)S^{p) 



(4.15) 



This equation is exact and it proves the Callan-Gross relation when we neglect the scaling violation (Q^-dependence) 

in the structure functions. 

Up until now, we can only trace the correspondence between our approximate equations (^|^) and ( [4.10D , and the 
LLA of the OPE-based calculations or the Lipatov's LL(l/x) approximation of the Regge calculus. The longitudinal 
structure function found via OPE has an extra small factor which is not compensated for by any big logarithm. 
Correspondence between the two approaches in the next orders is still unclear. In what follows we are going to use 
standard set of structure functions derived from the data in the leading logarithm approximations. Thus, we shall 
neglect all terms which can be explicitly reduced to F^. 



4-2. Renormalization of the evolution equations. 

Until now we have dealt only with objects which do not require renormalization. All these objects were tightly 
connected with observables. Corresponding Greenians and self-energies were imaginary and for this reason could not 
contain ultra-violet divergencies. So wc could safely use a "naive" form of Schwinger-Dyson equations, 

Gas = Gab + ^ Gar^rsGsb, Dab = Dab + ^ DarJ^rs^sb- (4.16) 

RS RS 

The divergent retarded and advanced self-energies were completely neglected, and retarded and advanced Green's 
functions were considered as the bare ones. Consequently, the counter terms of the Lagrangian still did not manifest 
themselves, and the renormalized coupling constant g,. still remains undefined. We must now fill this gap. 

After including the counter-terms we obtain the same equations, but with the self-energies modified by the quasi- 
local terms: 



Das = Dab+Y. ^^^[(^if H^s + ^iH^s) + (1 - Z:i){-lf5RsD^^]DsB , (4.17) 

RS 

Gab ^Gab+Y, Gar[Zif^rs + (1 - Z2){-1)'^5rsG^^]Gsb , (4.18) 



RS 



where 11' and 11" are the fermion and gluon loops respectively. In perturbative calculations the factors Ziaf should 
be split further as Zi^ip — 1 -I- (^i .if — 1) , with the second term assigned to the UV-renormalization of the vertex. 
The only changes from Eqs. ( 4.16| ) are due to additional diagonal terms in 11 and E. This is quite natural as the 



off-diagonal terms are im agina ry and if th ey were divergent we would have no remedy to cure the problem. As the 
matrix structure of Eqs. ( 4.18 ) and ( 4.18 ) is literally the same as that of Eqs. ( |4.16| ), we can rotate the 2x2 basis 



as usual Keeping in mind the light-cone dominance, we may rewrite equations for T-ordered and T^-ordered 

correlators as 

Doo(p) - D*{p) + Dret{p)[{,ZiFli[, + ^lU'/J ± (1 - Z^)Dq\p)]D adv{p) , (4.19) 



Goo(p) = G*(p) 



Gret{p)[ZlF^ll{p) ± {I 

00 



Z2)G^\p)]Gadv{p) , 



(4.20) 



and (2.11) with the bare vertices and no more 



where nii_oo(p) and Sh^qoIp) should be calculated using Eqs. (2^ 
than one source. We omit the *-labelled terms here as we are interested in the UV-renormalization and low- a; effects. 
These t erms are effect ively cut off at very high momenta and do not affect the latest stages of evolution. 

Eqs. (4.19) - (4.20) are too approximate to give an explicit value of running coupling constant. The questions 
we want to consider are: (i) does renormalization of the sources 11 and E result in renormalization of the coupling 
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constant? and (ii) does the coupling constant attached to the vertex at some moment t require for its renormahzation 
any information except for the dynamics of the ladder at previous moments? 

Before turning to explicit calculations, let us address some qualitative issues. The new equations are not completely 
independent from those we have already studied. Indeed, the sums and the differences 



Di = Daa + Dii = DiQ + Dou 

Do = Dret — Dadv = DiQ — DqI, 



Hi = Hoo + Hn = -Uw - 

Hq = Ilret — ^adv = ^ HlO 



Hoi, 



(4.21) 



define imaginary parts of T- and T^-ordered, as well as retarded and advanced, correlators in an interdependent way. 
The real parts are not independent either. Indeed, 



2Ds — Dqo — Dii — Dret + Dadv — 2Re£'oo 
2X1^ = IIoo — IIii = Uret + ^adv — SRelloo — 



-2ReDn = 2ReDret = 2ReDadv, 
-2Renii = 2Ren^et = 2Renad„. 



(4.22) 



If we recall, in addition, that the causality principle connects real and imaginary parts of the retarded and advanced 
correlators by means of dispersion relations, then we realise that we have practically no choice in the regularization 



of the divergent real functions - we must follow the BHPZ scheme |19 



Let us begin by writing down the integral equation for the gluon polarization correlators in the leading approxima- 
tion: 



n(;o(p) = n(^''(p) 



jTr[j^Gooik ~ p)YGret{k)[Z,F^,iip) + (1 - Z2)G^\k)]Gadv{k) 



-1-9r 



_(£k_ 



VaT/ ^ - -mice ik p)v:f;, {-p, p-k,k) 



x{Dret{k)[{ZiFKl + Zin'l,){k) + (1 - Z3)D^\k)]Dadv{k)}pf] (4.23) 

The corresponding equation for n^^(p) is the complex anti-conjugate of Eq. ( 4.23| ), and can be obtained via replace- 
ment the T-ordered functions by the minus T^-ordered ones and vice versa. Similar integral equations may be written 

for the fermion sources. 

In complete agreement with ( 1.21 ) and ( 4.22 ) these equations have ladder structure with retarded behavior. To 
lowest order, gluon and fermion Green functions are given by 



±d^'•'{p) 
p^±iO '■ 



Goo(p) = 



± 



(4.24) 



and n^^Q^^lp) is the usual ultraviolet-divergent vacuum gluon polarization tensor. 

For the sake of simplicity, let us consider only the gluon sector in the leading approximation. Projecting Eq. (4.24) 
onto the normal modes, one obtains 



where we have denoted 



2(27r)4y {k-py 



p+ wl^{k) + {l-Z3)k^ 



(4.25) 



I{z) - 8[( V/2 + k')Pgg{z) ~{k- pf{z + l/z - 1/4)] , 

and the splitting kernel Pgg is the same as in the evolution equation for a gluo n sour ce IIo i. It is IR-regularized in 
the same way, otherwise we would obtain a contradiction with equations (4.21) and (4.22). Alternating T-ordered 
and r^-ordered correlators in the l adder rungs is crucial for the subsequent conclusions. 

The imaginary part of Eq. (4.25) is finite. It can be obtained as the sum of the two ladder equations for vof^ and 



The divergent real part of the equation in the leading approximation is an equation for Rewj^ = Rcw^ 



the 



real part of the retarded self-energy. To facilitate the physical analysis, we shall derive this equation in another way, 
starting with the explicit expression for the retarded self-energy. Utilizing the identities Ilret = IIoo+IIoi — — IIio— IIn 
and Iladii = IIoo + Hio = —Hoi — Hn, we easily obtain: 



d'^k 



{2n) 



-p,k)Diir^{k^p)v^^^;,{-p,p+k.~k)D{ ''''\k) + 



-k - p. k)DT '"'{k + p)F;/;,(-p,p + k, -k)DiJf\k)] 



(4.26) 
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These equations have very clear physical meaning. The propagator in the loop is retarded (advanced) , and guarantees 
the required time direction. It is affected by the surroundings less than other correlators. The correlator Di describes 
the density of states which develops in course of the evolution. Thus, even the light-cone propagation is not really 
free - emission introduces additional phase shifts which result in the assembly of special wave packet. 



The sum and the difference of the Eqs. (4.26), respectively, are the equations for the real (lis) and imaginary (no/2) 
parts of the retarded self-energy: 

nr(p) = ^5.'/ ^KTf\p^~k-p,k)Df-'^^{k+p)v,:^;,i~p,p + k,-k)D('f^^^ , (4.27) 



nr(p) 



d^k 



cc ,a/3 / 



(fc) 



(4.28) 



Using Eqs. (4.2), ( 4.16| ), and ( 4.19| ), and projecting Eq. (4.27) onto the transverse normal mode, we arrive at 



2(27r)3 



w4i.Arri. ^21-^,P+^ ^^(fc) + (l-^3)fc " ^ 
d k5[{k-p) ]X( — ) — 



27r ] '^^^ {k-pY^'-k+' [fc212 



p+^w}(fc) 

>. V 7,-1- . 



(4.29) 



where V denotes princ ipal value integration. Another way to derive this equation is to separate real and imaginary 
parts in the Eq. (4.25). This consistency is a consequence of the dispersion relation for Ilret which allows one to 
recover Rennet via the already known Imllret- It reassures us that we are considering the propagation of the gluon 
field in proper environment of the pre-coUision cascade. 

One may easily see th at th e imaginary part of the self-energy, considered given, defines the free term in the 
inhomogeneous equation (4.29) for the real pa rt of the gluon self-energy. The unusual feature of this equation (which 
is common for all ladder- type equations like ( 4.23|) ) is the counter-term which is detached from the vacuum part of 
the self-energy. The latter is divergent, and the corresponding counter-term is in the integrand of the equation. This 
immediately requires that the kernel should act on the counter term as a 5{k — p). Futhermore, the integral which 
contains the counterterm, by inspection, is proportional to the one-loop vacuum self-energy of a gluon. After the 
UV-rcnormalization, the latter typically behaves like p^logijp' / h?), where A is the IR cut-off mass. We find that in 
order to have the counter-term — p^(l — Z^^) in its legitimate place near u'^Q-)-^(p), the following relation should hold : 

a,^l/log(pVA'). 

We have estimated the contribution of the counterterm by first calculating the imaginary part of the retarded gluon 
self-energy, and then singling out the logarithmic terms in the dispersion integral for the real part. We obtain 



ar (Pt) 



47r 



/3olog(p?/A2) 



(4.30) 



Thus, the known behavior of the running coupling constant is recovered. Since we have used very rough approxima- 
tions, the exact equality of the coefficient deserves further study. In the renormalization group approach it is a direct 
consequen ce of the Slavnov- Taylor identity in the null-plane gauge: Zi ~ Z^- After this renormalization the integral 
equation (4.2£) for the real part of the gluon self-energy [wl{p)Y'^" takes its final shape: 



wHp) 



w, 



(0) 



lip) 



(2^)2 



d^k5[{k - pf]T{ 



p+ wm as{pf)N, 
k+' [fc2]2 



(27r)3 



.P^.^(fc) 

(fc-p)2 [fc2]2 



(4.31) 



where the superscript ^^ren" is omitted. 

Despite the remaining uncertainty caused by the approximation, it seems to be very important that the running 
coupling has appeared as a consequence of causal evolution. Only the processes which took place inside the past light 
cone of the local interaction contribute to the magnitude of the coupling in its vertex. This guarantees the proper 
balance between propagation and emission-absorbtion processes at the pre-collision stage. 

We still have freedom to chose Z^. It has not yet been used in the renormalization; the running coupling has 
appeared as a necessary condition for renormalization rather than as an explicit choice of some physical parameters 
at some given 4-momentum. The strategy behind this choice must be the same as in the old-fashioned on-mass-shell 
renormalization of the asymptotic state: the on-mass-shellness means that the field propagation is steady, despite 
background vacuum fluctuations. Selecting this kind of boundary condition, we cannot describe the dynamics which 
leads to the "undressing" of the quark as required by the resonant condition of deep inelastic scattering. Moreover, 
the imaginary part of the self-energy must equal zero at the renormalization point. 
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Thus, the object we now study is a field config urati on which has quite different properties that those of a free 
particle. This is evident, for example, from Eq. ( 4.14 ), which indicates that imaginary part of w^{p) is strongly 



peaked near — 0. This configuration is singled out by two requirements (boundary conditions): (i) at the end of 
its evolution, it produces an off-shell quark that can interact with the electron in a resonant way; (ii) this quark stays 
bare (on-shell) after the scattering. 

In e-p deep inelastic scattering the second condition does not look too realistic. The bare quark will immediately 
fragment into a hadronic jet. In p-p collisions, the electron is replaced by a gluon from the second proton, but the 
final quark state must still propagate in the physical vacuum and decays into hadrons as well. Only in AA-coUisions 
do we expect the creation of a volume of perturbative vacuum large enough to allow almost stable free propagation. 

The most important point is that the pre-coUision dynamics of the field fluctuations is the same in all three cases. 
This is guaranteed by the geometry of the high energy collision and the causality principle. However, one should keep 
in mind that the type of fluctuations studied by DIS are strictly selected by the trigger of the specific measurement. 
Imposing other triggers will select other types of fluctuations. It can not be ruled out a priori that the pieces W2 and 
a2 that were neglected in LLA may become more significant. 

To find an analogy between the pre-coUision dynamics of the proton constituents and physics of continuous media, 
we may try associate components of gluon self-energy with electromagnetic susceptibility. If its imaginary part is 
infinite then we have an ideally conducting medium. (Remember that p~ is a "frequency" corresponding to the 
"time" .) A significant growth of the imaginary part of the self-energy at the "resonance" p~ = should lead to 
anomalous dispersion - the real part must drop to zero : 

Reu;f(p2,p+,p-^0)-(l-Z3)(-p?)=0 . (4.32) 

(We can suggest the formal mathematical argument: once Imllret ~ ^{p~), then from dispersion relation Kellret 
V{l/p^) which, though in a singular manner, is equal to zero at p^ — 0.) In the region of anomalous dispersion, 
the phase and the group velocities should have the opposite signs. This reveals one more unusual feature of the 
"undressing" process: the phases of the fields participating in the assembly of the wave packet representing an 
interacting quark (or gluon) travel in the "normal" time-direction (from the past to the future). These fields leave 
"holes" in the sea. Only the act of measurement (scattering) transforms the creation of these "holes" into the process 
of multiple emission. 

5. Distributions of quarks and gluons in leading order. 

In this and following Sections we present results of an explicit calculation of the single-particle distribution of light 
quarks and gluons produced at the earliest stage of an AA-collision. In this Section we calculate cross-sections to the 
lowest order. It is general practice to associate the corresponding processes with the excitation of the sea quarks and 
gluons. There are three processes of this type. Their graphs are shown in Fig.|l| In Section 6 we will derive a complete 
set of equations for the first order processes, and discuss how one avoids difficulties which accompany calculations 
based on the factorization theorem 



5.1. Production of light quarks to leading order. 



We return to the initial formulae (2_^) and (2_^), and rewrite the former in the momentum representation 

dN, _ TT[i i>^^M] 



(27r)32p" 



(5.1) 



The lowest order of our theory assumes: (i) the exact vertex operator must be replaced by the bare one, which leads 
to 



P 







9 



dpd'^x 2(27r) 



Tr[i,r7^Goi(p + fc)tVD?g,,^(fc)] 



(5.2) 



and (ii) in Gqi and Dio possible contribution of the out-states of quark and gluon in the continuum are excluded. 
Thes e co ntributions are described by the higher orders of the perturbation theory. Thus, wc naturally arrive at 
Eq. (|~ 



dN, 



dpd^x 2(27r)3 



^^S{k + g -p){Tr[#^''G(f' (9)t''7^D(t),^;(fc)] + (A) ^ (B)} 



(5.3) 
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were the quark and gluon correlators must be taken in the foUowing form: 



*^01 - ^-'Ol - ^ ret ^01 *^adv^ ^01 "^01 ' ^ ret ^^01 ^adv^ [•J-^,!^), [bA) 

with Gq\ and DqI representing the quark and gluon distributions at some arbitrary scale Qq. For computations, 
we shall use a standard CTEQ parameterization of the nucleon structure functions [gO| q{x,q^) and G{x,k^). These 
were obtained by fitting the data with the solutions of the GLAP evolution equations. For nuclei, we shall also use 
the semi-empirical formula fit to nuclear shadowing data. 

We proceed in the laboratory frame, which is the infinite momentum frame for both proton A and proton B. The 
directions of their light cone propagation are fixed by the two null-plane vectors and n^: 

n^ = (l,Ot,-l), 71^^(1,0^,1), n\=nl=Q. (5.5) 

They define the light-cone components of the Lorcntz vectors: 

uacl — = = oP + a^, 71^6 = 6^ = 6+ = 6° — 6^. 

We split the total cross-section into the three parts: 

<) =<)(v,n) + 4")(g,i]) + a(")(n,i]), (5.6) 

The term crg*^''(tj, t/) is naturally absent, as energy-momentum conservation prevents fusion of the two on-mass-shell 
particles into one on-mass-shell quark. 

Unlike the case of DIS, both invariants from quark (ai and <72) and gluon (wi and 1x12) self-energies contribute to 
the cross-section of quark production. In order not to exceed the accuracy of the leading logarithmic approximation 
of the GLAP-evolution of structure functions, we shall omit <T2 and W2- They are of the next order by a formal count 
of the a^-powers, and are not under direct control of DIS data. 

Any term coming from the nucleus A carries a 6{k~), and any term coming from the nuclei B carries a S{q'^). This 
drastically simplifies the calculation. The first two terms from Eq. ( |5.6| ) are calculated explicitly, yielding 



— -rn — = ^ — [G{—/^^Qm{—r^Pt)[^ + J + (y^-y)] , (5-8) 

dpfdy 3s y/s y/s sxaxb 

where y is the longitudinal rapidity of the final state quark, xa,b = Pt&^^ I \fs, and the substitute term (y —y) 
accounts for the symmetry {A ^ B). We have also defined G'{x,p1) = dG{x,p^/dpf and q'{x,pl) = dq[x , pi) / dpi . 
The main contribution to the cross-section comes from the third term in Eq. (5.6): 

d4''^(E,n) _ Svrao [ dk1dq1{k1 + q1){l + ql/sxAXs) [g^( P*^~^ q^)G' {^^^^ k"^) I (y : y)] (5 9) 



dpidy ispl J ^[{kt + qtY-plM-ih-qtY] Ws 

where all terms which are at least as small as the longitudinal structure function Fl of the DIS were neglected. The 
square root in the denominator comes from the angular integration over the orientations of the transverse momenta: 

/dWq*<5(^Hk, + q,-p,)F(fc*,gO= /^^^^^^' (5.10) 

J J 4:S[kt,qt,Pt) 

where S{kt, qt,Pt) is the area of the triangle with sides fcj, qt and pt, 



^S{kuquPt) = ^[{kt + qt?~p1][p1 ~ {kt - qtV] , (5.11) 
and the integration domain is restricted by the triangle inequalities, 

kt+qt>Pt, \kt-qt\<Pt ■ (5.12) 



19 



5. 2. Production of gluons to leading order. 



For the case of gluon production, we start from Eq. (2.9), and, as for Eq. ( |5.l| ), rewrite it in the momentum 
representation : 



dpd'^x 



df,,.{p,u)[-inQ'^'^''{p)] 



(27r)32pO 



where the projector dfj_^{p, u) 



pt'u" + u'^p" pf'p" 



(pu) 



{pu)'' 



, d'^''(p,wK = 0, p" 



(5.13) 



(5.14) 



is a sum over the two physical gluon polarizations in laboratory frame with time axis along the four-vector u^^ 
(1, 0, 0, 0). For this gluon we choose the time-like axial gauge m''_B" = 0, = 1. 



To the same approximation as in Eq. (5.2), we may write the gluon polarization tensor (temporarily omitting the 
fermion contribution) as 



n[;r (P) = -^9l{ I ^J^S{k + q - p){VZ^ik + ~q, -fc)D(f Ji?(fc)n^,^i! i-k q, q, fc)Dlf ^/^(q)] + {A 



B)} (5.15) 



The gluon distribution of each nucleus consists of the two familiar terms given by Eq. (5.4). Again, we shall keep only 
the leading terms in the gluon sources which are under the control of the data : 



[d{k,nA)Il'^''Hk)d{k,nA)]Z' 



add' 



-d'"'{k,nA)wt{k) 



(5.16) 



A similar expression may be written down for nucleus B. Calculation of the trace over the abundant color and vector 
indices results in an overall factor of 24 x 32 x T{k, q). The trace T(fc, q) is given below. As in Eq. (5.6) we have two 
types of contribution to the total cross-section : 



4°) = a(°)(^,n) + a(")(n,n), 



The first term comes from the interaction of the "source" with the field from the "initial" distribution 



dptdy 



y 



t)G\^,pi) 



y 



1 



ple^y 



s{xA + xbY 



+ {V'^ -2/)] 



The second term comes from the interaction of the two sources : 



d<jf\n, n) ^ 24^Qo r dk'^dq^T{kt,qt 

dPtdy spj J 4S{kt,qt,pt) 



G'{?^,k^)G'{ 



Pte^ 



The trace T{k, q) is conveniently written in the following way 



T{k,q)^k^ +q^ + 



2^2 



(kf+qf) 

S{XA + XbY 



qt 



kWt 



H ^ 

Sx\ S^XAXB 



'2kWt 



s{xA + xbY 



ql 



kl 



(5.17) 



(5.18) 



(5.19) 



(5.20) 



The integrand of Eq. ( |5.19| ) is symmetric with respect to interchange of fcf and qt. H ence the contribution from the 
interchange term (A <-> B) is the same, and factor 2 is already included in Eq. ( 5.19 ). 

We close this Section by discussing two minor problems. Th e firs t is an apparent infrared divergence of total 
cross-section due to the kinematic factor p~^^ in Eqs. (|5.9| ) and ( 5.19 ). However, one can easily see that this factor 
disappears when the integration variables kt and qt are rescaled by pt. After such a rescaling, any possibly remaining 
issues concerning low pt behavior are connected with the sources (or structure functions) which are controlled by the 
data. 

The second problem arises from the positive powers of transverse momenta of the incoming quark and gluon fields. 
This can be seen from Eqs. ( |5.9| ), and (5.1£) and (5.20), and appear in the same form in the next perturbation orders. 
The formal solution of the problem is very simple: to be consistent with the condition under which the sources (or 
structure functions) are defined, viz.., s oo, we must drop all terms containing ratios like qt / sxaxb and k ? /sx axb, 
despite the external kinematic condition sxaxb ^ Pt- If this is not done then the integrals (5.£) and (5.19) will 



20 



diverge, as the LLA structure functions do not provide sufficient cut-off at high kt and qt- In the context of the 
existing theory, these ratios are parametrically small and may not be switched on until the high-order twists are taken 
into account. 

Nevertheless, the growth of the matrix element at high transverse momenta is quite physical, and may be interesting 
at very small x. The wave packets with small x are smooth and extended in the longitudinal direction. If we require 
that large kt and qt add to form small pt , then the initial geometry of momenta is almost coUinear and the incoming 
fields effectively overlap. We thus encounter a type of coUinear singularity which may be shielded only by appropriate 
behavior of structure functions. The latter is naturally provided by the BFKL equation (See Ref. |10| and Appendix 
1.) We emphasize that this process may be taken seriously only in the domain of very low x, where only very 
preliminary data exist at the present time. 

6. First order corrections to the quark distribution. 

In the previous Section we have calculated the lowest (Born) contribution (of a^'Oi'dcr) to quark and gluon pro- 
duction corresponding to the "excitation" of a single quark or gluon from the sea. Now we intend to show that the 
next order terms of our kinetic perturbation expansion describe the more usual "creation" processes. They do it 
in a very special way such that no infrared divergencies appear in the calculation of the total cross-section. This 
may be the most important result of the present calculation. Furthermore, this requires no special proof: because 
the perturbation series for the probabilities was initially resummed in the new expansion, it does not generate any 
terms affected by the initial state coUinear singularities. These are absorbed into definition of the sources (structure 
functions), and thus solving the problem of divergence in an experimental way. 

This is in striking contrast to all known calculations based on the factorization theorem. Such calculations always 
start from a master formula, such as 

yaiAB-^pX) ^^^ /'d.. f\xy-^^p^FM^^,Q^)FM-,,Q^) , (6-1) 

P a,b ceX P 

which heavily relies on the parton model, and implicitly impart the status of an observable at least to one additional 
final state particle. The 2-^2 processes have the lowest order in this approach, and the Q^-dependence of the 
structure functions Fjj defines the so-called factorization scale, rather than reflecting its full QCD-evolution. As was 
already mentioned in Section 4, the reason structure functions are treated mistrustfully in this approach is connected 
with the out-of- light-cone behavior of Feynman's propagators for the incoming partons. 

Quite naturally, apart from all other amplitudes, the 2 — )■ 2 cross-section includes those corresponding to the 
emission of a second particle c G X from the initial state. The squared moduli of these partial amplitudes duplicate 
those already included in the definition of the structure functions at lower factorization scales. 

The new approach does not distinguish any states absorbed into the set X. At the same time, it restores the proper 
status of the retarded propagation for the incoming fields. As a result, it excludes a priori any double counting of 
processes. 

The first order corrections to the quark or gluon production can be divided into two major categories: the self- 
energy-like and the vertex-like. These names reflect only the topology of the new diagrams. 



6.1. Self- energy-like terms. 

The sequence of the cts-order diagr ams in the sclf-energy-type term emerges from the possibility that cither the 
quark of gluon field correlator in Eq. (5.2) represents a field in the continuum of out-states. 

J AT „2 r rt^k 

P'^dk = 2(2^ J (^{Tr[i,r7^Go# (p -f k)t'>r-D\l,^{k)] + Tr[^t^rGoi{p + k)t\'^ D*%{k)]} , (6.2) 

This immediately means that the remaining exact correlators Gqi and Dig carry information about both nucleus A 
and nucleus B. In other words, these fields were created during the course of the collision between the two nuclei. The 
only terms of the subsequent expansion of the quark and gluon correlator that survive in this case are 

Dio = -D^nioD*,, and Gqi - -G#,SoiG#, (6.3) 

The superscript indicates that we consider propagation after collision of the nuclei, and in case we need radiative 
corrections to this propagation, we must consider them against the background of the distribution of quarks and 
gluons created in the collision itself. 
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The intensities 11 and S of the field sources created by the two nuclei, were already calculated in the previous Sec tion , 
but now these fields are off-mass-shell, and the terms with products of two ^-la bell ed correlators from equations (5.4) 
should be added. These are just the lowest order terms of the master formula (6.1) of the QCD parton model, with a 
fixed factorization scale Qq, and correspond to processes of the type 2^2, where parameters of the second emitted 
particle are completely integrated over. Thus the complete formula will be 

+Tr[i,tV[G#,(fc)I]^i^(A:)G#„(fc)]tV^fo'^^(fc-p)]} , (6.4) 
where the following chain of substitutions is supposed : 

^ofiP) - *.90 / (04'5(fc + 9 - P){^V[Go^l*(fc) - Gr^m^Olimad.im'j^ X 

iD^O*iq) - T)retiq)lltoiq)'Dadviq))ba,.f. + iA^ B)} (6.5) 
for a "collective" quark source, and 

V [Df„*(fc) - D.,,(fc)nfo(fc)D„,„(fc)] VV[Dt,*{q) - Bret{q)Uto{q)T)aM] + 

+TYtj[G^{{k) - Gret{moiik)Gadv{k)MGf^^-{q) - Greti^^^ , (6.6) 

for the "collective" gluon source. Here V-s stand for the 3-giuon vertex with all arguments dropped. 

Three of the six graphs of this type are given in Fig.^. The first corresponds to the s-channel part of the Compton 
process. The other two are of an annihilation type. They can all be represented as squared moduli of the corresponding 
amplitudes. The t-channel partners of these diagrams are not generated by this perturbative expansion. If they did, 
they would duplicate processes already included in the definition of the sources. The three omitted graphs correspond 
to the interchange (A <-> i?). 



6.2. First order terms from the "vertex-like" corrections. 



The general expressions for the one-particle quark and gluon distributions (Eqs. ( |2.8| ) and (2.11)) contain dressed 
quark-gluon and three-gluon vertices 

y',^lHSB(^:V,y) = (-i)«+^+-1P:^|MIl!. (6.7) 

Because of their additional matrix structure they contain some unusual elements which should be determined now. 
It can be done easily by using the formal solution of the matrix Schwinger-Dyson equations : 

[G-^]ab = [G-%B - ^AB, [T^-'Ub = [D^^Ub ~ TIab. (6.8) 

Functional derivatives of the bare Greenians give us the bare vertices, and those of self-energies give the corrections. 
Using the one-loop formulae for the self-energies it is straightforward to find the first order corrections to the vertices : 

^'^r^B,5(e,y;'7) = ~^9'r{~lf'+'{l't'^Gns{i,0l't''GsBiC,yh''t'B';^^''{y,0 + 

+7''t-G«,B(f, yh'^t''B'j^-f{y, O^.^l (C)D^T'(C, 0} • (6-9) 



Eq. (6.E) contains an additional rule: besides the expected contour indices, the correction ^^•'F^^^ g(^, y; rj) to the bare 
vertex with the contour index B acquires an additional factor (— 1)^^"'^. This rule works in the same way for the 
three-gluon vertex which is too cumbersome to be written down separately. 

The sum over the contour indices naturally divides into two groups: those with R — S and R ^ S. We delay 
discussion of the first group which is responsible for those types of radiative corrections which may lead to vertex 
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form-factors. The second group, which we will discuss below, describes contributions of the higher-order real processes, 
like emission of a second uncontrolled jet. 

The one- loop vertex- type correction to single light quark production contains four terms, two from the first term 



in Eq. (3.£) and two from the second one 



where 



Xi = t"7^Goi(fc)tVGio(fc + q - p)t VGoi(g)tS''Dgg^^^(p - k)Dil,Jp - g); 
I2 = ^^7^600^ - fc)t VGoi(p ~k- q)tVGii(p - <z)tVDgJ,^p(fc)DgJ,,,(g); 

X3 = f^T^Goob - fc)t=7"Goi(<z)tS^Dg-;^(fc)D-;,,(fc + q - p^l^^ip - q)Vf,U-k, k + q-p^p-q); 

X4 = t"7^Goi(fc)t^7-Gnb - q)t'Y'D^oCpiP - k^ixAP - k - q)B'^oiM'lK'lUk -p,k + q-p,q) . (6.11) 

These four terms are depicted in Fig.^. There are three cut lines in each diagram, and they can be "distributed" in 
six different ways between the two colliding nuclei and the additional out-state excited in the continuum. After that 
we can convert every term of this expansion into a product of two scattering amplitudes. This is done in Fig.^ (up 
to a trivial interchange A <^ B). It is clearly seen that amongst this group of graphs there are none which would 
represent the squared moduli of any amplitude, but all allowed interference terms between all the processes which 
produce a quark and something else are included. 

Recalling the previous discussion of the self-energy-type terms, we see that our perturbative expansion docs not 
generate any diagrams which would repeat those present in the definition of the sources (via their ladder expansion) . 
These missing patterns are not IR-safe, and were regularized and renormalized in the course of their definition via 
the DIS cross-section. 

We shall now show that no further infrared problems appear. We demonstrate this using the definite subprocess 
of a detected quark with momentum p accompanied by an uncontrolled anti-quark in the out state. They were both 
created in a collision of two giuons, gA9B ^ 99 • 



6.3. Infrared safety in the Us- order. 



Infrared finiteness of the self-energy-type diagrams of Fig|| is intuitively understandable, since the intermediate 
s-channel gluon or quark carry a large time-like momentum. We may expect IR-problems only in the vertex-type 
diagrams, like (VI) of Fig.^, where the intermediate fermion is in the t-channel. The corresponding distribution of a 
single quark is given by the expression : 



dpd-^x 2(27r)3 



Dii>ik) - B^->{k)U\->ik)r>^-i{k) D\r{q) - Bl^:>{q)Il\-'{q)B\2iq) 



d^kd'^q 
(27r)8 



Tr[#'^7^Goo(p - fc)tVGoi(p - fc - qX^GniP - 9)tV] 



JB)* 



bd 



(6.12) 



where retarded and advanced functions carry a superscript indicating which nucleus was the source of field. The 
natural variables for subsequent calculations are quark rapidity and momentum fractions defined via 



P 



Pte 



±y 



k+ 



^SXA, q 



The most troublesome element in all following calculations is the trace over spinor and vector indices, 

T - Tr[i>7^(i.- h~ h)Y{t^ A)Y]dp,{k, n^jd^q, ub) . (6.13) 

As in the lowest order, we can single out different types of terms contributing to the cross-section in the first order: 



aW = 41) (g, g) + (fj, n) aW (n, n) 



(6.14) 



The first term corresponds to the first nonvanishing order of the master formula of the parton model, which factorizes 
the "hard" QCD cross-section and structure functions evaluated at some (sufficiently high) scale Qg. Two mass-shell 
delta-functions make the calculations relatively simple, and the result reads as follows: 
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d<7q^\G,G) _ TraQ dxAdxB 



dpjdy 125^ Jo XAXB 



SXAXB 



G{xA,Qt)G{xA,Qi) X 



x5[xAXB - "^{xAe-y + XBey)]e{xA - ^ey)0{xB - ^e-y) . (6.15) 



We see that this term remains finite when pt and is strongly suppressed, both by the second power of s in the 
denominator and the smallness of a^. 

To facihtate the following analysis, let us trace how the structure emerg es. W e may expect an infrared divergence 



from the poles of the product of the two i-channel propagators in Eq. ( 6.12 ), Gqq{p — k)Gii{p — q). Since they 



arc unshielded by finite masses or virtualities, the factor (p^sxaxb) appears. One power of s comes fro m the 



definition of cross-section, and one more from the final-state phase space. The combined spinor- vector trace in (6.14) 
gives a factor p^s. As a result, no large logarithm which could partially compensate for the smallness of the coupling 
constant appears in the total cross-section. Eventually, we may expect only a weak scale dependence of the total 
cross-section. 

The "box diagram", which is an unavoidable partner of crq^\Q,Q) in the approach based on the factorization 
theorem, is IR-divergent, but it has no analog in our perturbation expansion. Nevertheless, we considered it useful to 
calculate it in Appendix 2, in order to compare its structure with that emerging from the present calculations. 

The mixed term cTg^^ (^: n) , and the term cTg^-' (11, H) (contributed to by two "sources" ) are more complicated because 
one or both of the incoming fields interact with their sources and hence are off the mass shell. This means that at 
least one of the mass-shell delta functions is no longer present, and traces become unwieldy. We will not present their 
explicit form here, as these terms are not expected to be large. Indeed, after the box diagram has been extracted, 
one retains only interference terms. These are usually small, unless they are affected by a singular infrared behavior. 
A sufficient qualitative analysis of this behavior can be done without explicit calculations. It is enough to notice that 



the trace T(p,k,q) in the numerator of the integrand of Eq. (6.12) is, in general, a polynomial of fourth order with 
respect to pt/y/s. When the momenta of both structure functions are put on mass shell, the only powers that survive 
are pf/s and pf/s"^. This leads to an immediate cancellation of the two unshielded infrared poles of the propagators. 

In the integral for aq^^ (11, 11) both poles are shielded by gluon virtualities, and infrared divergence can not appear at 
all. In the integral for aq^^Q^H) only one pole is shielded, while the second produces an unpleasant p^^ behavior. 
However, this does not lead to an infrared divergence of the cross-section, since the same power p^^ is implicitly 
present in the differential dp^ on the l.h.s. 

In Appendix 2. we consider in full the s-channel production of light quarks from the process gg — > qq. The 
mathematical details behind the above qualitative analysis can be found there. Here we shall discuss only the main 
physical issues. 

The higher powers of pt in the trace (|6.13 ) lead to an increase in the differential cross-section at high pt over the 



lowest order result. This is in compliance with the observation that the first order terms bring more than a simple 
quantitative correction to the lowest order. It is precisely the emission of two back-to-back jets which makes possible 
the existence of a high-pt particle in the final state. 



7. Conclusion. 



We have considered new principles to compute the distributions of quarks and gluons created in the first hard 
interaction of the two heavy ions at high energies. We essentially employed an initial resummation of the perturbation 
series for the probabilities It allowed us to describe two different high energy processes, viz., e-p scattering and 
nuclear interactions, in the same terms, as two versions of the same phenomenon - deeply inelastic scattering of 
composite systems. 

It is shown that these calculations can be performed without reference to parton phenomenology. We have in- 
troduced the concept of a source as the main subject of QCD evolution, and have shown that the equations which 
describe the dynamics of the sources are independent of the type of high-energy process, and independent of the 
particular choice of the final interaction. 

The additional benefit of the new approach is that it explicitly displays the causal structure of the QCD evolution 
equations, and their physical meaning as the spectral analysis of the composite system as performed by the interaction 
which results in the bare quark or gluon production. The evolution equations for the sources allow for a smooth 
transition between the regimes described by the GLAP and BFKL equations. The by-product of this study is a new 
form of the fusion term in the GLR-type equation, which might lead to the stronger low x saturation of the sources 
than any terms considered previously. 

One of the most important results of this study is the new type of perturbation expansion, which, unlike for the 
factorization technique, does not lead to double counting of processes. The diagrams already included in the definition 
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of the sources, and controlled in aggregate by the DIS data, do not appear again in the higher orders of the new 
perturbative expansion. The diagrams which do appear are free from initial state infrared (coUinear) singularities and 
do not require artificial cut-offs. The leading parts of these diagrams do not depend upon the initial factorization scale 
either. The price one pays for the efficiency of these calculations, is that one requires the full x and dependence 
of the sources (or structure functions) extracted from the data. 

We are now in a position to calculate the single-particle distributions of the quarks and giuons after the first 0.1 fm 
of the heavy ions collision. 
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Appendix 1. The full form of evolution equations. 



In the main body of the paper, we explicitly studied only the part of evolution equations which eventually results 
in the GLAP equations . T he latt er a re also known as the leading logarithmic approximation (LLA). 

Corrections to Eqs. (|4.3| ) and (4^) are of various origin. Unfortunately, we cannot count these corrections in the 
traditional way, which relies on the firm hierarchy of twists in OPE-based calculations. For the moment, by "next 
approximation" we shall mean some kind of structural expansion based on the complexity of the processes taken into 
consideration. Surprisingly, it does not lead far away from the commonly used scheme. 

In what follows, we study two types of c orrect ions . First, we consider the terms missing in the analysis of the 
simplest (by their structure) equations, (4.3) and (4.4). In the next two subsections we study the low- a; region which 
results in BFKL equation [|l0| as the limit of new evolution equations. More complicated corrections lead to an 
equation resembling the GLR equation pll,^2|, but with some significant differences. 



Al.l. Trivial corrections. 



The simplest corrections arise from Eqs. (O) and (4^), as a residue of the original spinor and tensor form after 
LLA-terms have been extracted. To make them more visible, let us split <Ji{p) and Wi{p) into a leading ( a'l and w'l) 
part, corresponding to LLA, and subleading parts (ct", w'l,...) : 



<Jiip) = <ip) + <{P) + <{P), w^{p) = w[ip) + w'lip) + wTip). 



(Al.l) 



The mathematical steps just follow those of the perturbation calculations described in Section 3. Now it is only a 
lengthy exercise to obtain the resulting equations. Spinor components of the quark source may be expanded as follows: 



a[{p) = / d^kA 



-p^k^ 



p- 



k+ai{k) 



fc+'5f(fc)5f(/fc) 



p+ fc+Wi(fc) 



fc+'Wf (fc)Wf (fc) 



(A1.2) 



f^iip) = J d'^kAkpk'^ 



p^ ^ fc+(Ji(fc) 



k+^ S^ik)Sf{k) 



k^wi (fc) 



fc+'Wf (fc)Wi^(fc) 



(A1.3) 



a'I'ip) = / d^kAkpP^ 



' 0-2 (fc) , p+ k^W2{k) 

^S^{k)S^{k)^ k+'w^ik)Wi{k) 



(A1.4) 



where we have denoted 



Akp — 



Cp 



k+(7i{k) 



(27r)3fc- 



:S+[{k-pf] = 



2gfeik"~p°: 



k^wi{k) 
k+'W^{k)Wf{k) 



(A1.5) 



S[{p+-k+){p- 



k-)^{pt-ktf] 
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The meaning of this decomposition becomes clear if we muhiply it by the kinematic fact or p"*" , a nd sa ndwich it 
between retarded and advanced propagators (see Eq. ( 4.13| )). The joint left hand side of Eqs. ( Al 2 )- ( Al.4 ) becomes 
the derivative with respect to of the quark structure function of DIS. The r.h.s. of Eq. ( Al.2) then leads to the 
GLAP pa rt of the evolution. The factor in it is responsible for the logarithmic behavior. Now it is easy to see that 
r.h.s. of (Al.3) will result in a term with l/Q^ behavior. It simulates the next twist contribution, even though the 



second twist is not in cluded explicitly in the density matrix as an effect of next order correlations. 

The r.h.s. of (Al.4) represents a secon dary i nfluen ce of the longitudinal components of the spinor and gluon sources. 
This components are defined by (A1.5) and (Al.8), and are of next order by a formal count of the a^-powers. a-i 
directly contributes to the longitudinal structure function given by Eq. (4.15). It also appears in the Born term of 
the quark excitation process, but we leave it aside for now since it is poorly controlled by the data. 

All of the above comments apply equally to the components of the gluon source, which can be decomposed following 
the same principle: 



d*kA 



kp 



p+ fc+cri(fc) 



fc+'5f (fc)5f (/c) 



p+ fc"*"i/7l(fc) 



fc+'Wf (fc)Wj4(/c) 



(A1.6) 



d'^kAkpk 



p+ k+aiik) 
'''^k+'s^{k)St{k) 



p+ k+Wl{k) 

''^k+^W^{k)W^{k) 



(A1.7) 



w'I'ip) 



d'^kAkpk+ 



<y2{k) 



k+'smsHk) 



ON (?— ~ ^2 k'^Mk) 



(Al. 



W2{p) 



= Jd^kA^.'^i-SCFUfil-^)-^ 



k+cri{k) fc+ _ p+ 2 k+wi{k) 

{k)S^{k) 'p+^ 2k+' W{^ik)Wf{ky 



(A1.9) 



A preliminary examination of the additional terms in the evolution equations reveals that they have the same 
singular infrared behavior as GLAP equations, and should be regularized and renormalized. The way to do this is 
not yet clear, since conservation of momentum seemingly fails to control all necessary subtractions. A complete study 
of these equations is a separate subject. In this paper, we have considered the explicit form of the DIS structure 
functions as granted. 



Al.2. The BFKL equation. 



As has been shown in Sec. 4, the new evolution equations for the sources, even in their reduced form ( |4.3| ) and ([4.4|), 
are practically the equations for the derivatives of the structure functions, rather than for the DIS structure functions 
themselves. The former were first introduced by Lipatov as the "unintegrated structure functions," and they relate 
to the latter in the following way 



2^ _ 2d x G{x,p1) 
dpi 



Thus, up to an insignificant normalization factor, we may write: 



W^{p)W^{p) 



(ALIO) 



(ALU) 



In the limit of low x, the function f{x,p1) was proven to obey the so-called BFKL equation, 



df{x,pl] 
dx 



3«s(Pt) 2 

.10 



d^ 

k? 



f{x,k?)~f{x,pl) , f{x,pl) 



\pl~k?\ 



P\ 



(A1.12) 



This equation was originally obtained by considering the amplitude of the process 2 — > 2 + n and summing the leading 
log[\lx) terms | |To| ]. 

Evolution equations like (^^) were derived immediately as integral equations. At high and not too low x, 
they allow for simplifications which lead to the GLAP equations. This means that the necessary resummation of the 
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leading logQ^ was done from the very beginning. Next, our intention is to determine what simpHfications should be 
done to obtain the BFKL equation. These simplifications are effectively equivalent to a reduction in the number of 
diagrams already accoun ted for in the integral equa tion. 

In the notation of Eq.( AI.IO ), our equation ( [4.10|) for the transverse gluon source reads 



Sip )f{x,pi] 



TT'^Pt I — 



y 



6[{p+ - k+)p- - (p, - kO^] 



P,,(-)/(y,fc?) (A1.13) 



We have retained all terms that contribute to the LLA. We notice that at real momenta the denominator 
yVi {p)yVf'{p) is strictly positive both before and after any integration without additional weight. Thus the full 
integral is not singular. The next steps are as follows: 

Step 1. For x ^ 1, approximate the splitting kernel as Pgg{x/y) ~ (y/x). 

Step 2. Integrate both sides of the equation over p~ , considering the retarded and advanced propagators in the 
integrand on the r.h.s as bare. In this way the singular infrared behavior of the integrand is unshielded. Integrate 
over the azimuth angle. 

Step 3. Expand the propagators on the l.h.s. up to the first order in radiative corrections. Retain only the vacuum 
correlator and move it to the r.h.s. in order to shield the singularity that resulted from the approximation. 
Step 4- Differentiate both sides of the equation wi th resp ect to x. 

This procedure will result in the BFKL equation ( Al.12 ) (up to the insignificant last term in the integrand). The 
asymptotic behavior of its solution is well known. Its exponential part, 



f{x,Pt) ~ exp 



2A" ln(xo/a;) 



(A1.14) 



provides a decrease of f{x,pf) at high that is faster than any negative power of p^. 

We postpone any discussion of the accuracy of the above approximation, and satisfy ourselves with the most 
important fact that the evolution equations contain both GLAP and BFKL regimes of evolution as the limits. Thus, 
we may hope to describe a smooth transition between them. From a pragmatic point of view, the above asymptotic 
behavior guarantees the convergence of the integrals that appear in a calculation of the cross-sections of quark and 
gluon production (see Sections 2 and 6). 



A 1.3. Gluon shadowing. 



Mo re c omple x co rrections to the GLAP evolution correspond to the replaceme nt o f the ^-labelled correlators in 
Eqs. (L3) and (L4) by correlators with sources, viz., the second terms in Eqs. ( |4.2| ) and (L2). This means that 
instead of the final state correlators describing the emission, we include the initial state correlators. This replacement 
accounts for possible fusion of the partons. Fusion is expected to be most important for gluons at low x. So, only the 
purely gluonic component will be considered here. 

Terms responsible for fusion of two gluon fields have the form 



/ ^^<5(fc + g-p)lC7(fc + g,-g,-fc)x 



X [Dret{k)Tlio{k)Dadv{k)]ff. <,^/,(-fc - kk) [Dret{q)^w{q)Dadv{q)] 



pX 



(A1.15) 



This equation is identical to Eq. ( 5.15| ) describing gluon fusion in a nuclear collision, except that here both gluons 
are taken from the same nucleus. We shall use the standard approximation, i.e., bare tree propagators and bare 
vertices. The longitudinal gluon function W2 will be neglected also. Routine calculations similar to those performed 
in Section 5 then yield: 



AfusG'{x,p1) 



37rao (27r)3 ^ 
' 2pt ttR^ 



XlX2 



XiX2 



X^ 



dxi r dx2 , . 

/ d[Xl+X2-X) 

Xi Jq X2 



dkfdq^ 



AS{kt,qt,Pt) 



1 

Xl 



X2 



pI 



G'{x,,k?)G'{x2,q^), 



(A1.16) 



where G'{x,p1) = dG' (x , pf) / dpi and Slkt, qt,Pt) is the area of a triangle with the sides kt, qt, and pt. The initial nor- 



malization factor {ViabP^) ^, which was convenient for calculation of cross-sections, has been replaced in Eq. (Al.K;) 
by (7ri?2)~^, which corresponds to a normahzation per unit transverse area of a nucleus with radius R. 
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Eq. ( A1.16| ), by its structure, is very similar to the well known Gribov-Levin-Ryskin (GLR) equation ||T^, and one 
derived later by Mueller and Qui jlj]. It clearly reveals the same tendency to saturate the rate of the field source 
QCD-e volutio n at low x. However, it seems to have several differences. The most significant is that the power of as 
in Eq. ( Al.f 6 ) is less than in Refs. [0 and jl^. The formal reason is simpler form of the vert ex of t he " three-ladder 
interaction" that is prescribed for us by the general structure of the evolution equations (4.3) and (4.4). The other 
differences are of dynamic origin and will be discussed elsewhere. 



Appendix 2. Some estimates of the first order terms. 



A2.1. The "box" diagram. 



Tw o box-type diagrams appear in the calculation of the quark production cross-section, if we use the master formula 
( |6.l| ). They are depicted in Fig.^ and the corresponding analytic formula is : 



2pO ^ 



dpd'^x (27r)3 J (27r) 



d^kd'^q 



[Tr #^''G„t(p - g)tS-G# (p - k - qXY 



Routine calculations result in the following expression for the differential cross-section : 
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{box) 
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dptdy 
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xaxb Vs Vs 



p^(3 - xaxb) 
sxaxb 



Apt 



{sxaxbY 



G{xA, Ql)G{xA, QI)S[xaxb 



Pt 



{xac y + XBey)]. 



(A2.1) 



(A2.2) 



The factor p^^ appears in the same way as in Eq. ( ^.9D for the Born term. H oweve r, previously it could be effectively 
absorbed into the structure functions, which is not the case now. Indeed, Eq. (A2.2) has no additional integration over 
momenta which could be rescaled by pt- We expect at least a logarithmic divergence of total cross-section as a result. 
This divergence may be strengthened by the low- a; behavior of structure functions because of the pt dependence of 
the low limits of integration over xa and xb in Eq. (A2.2). 



A2.2. Inclusive production of quarks in s-channel. 

The analytic expression for the diagram (S2) of Figj3[ corresponding to the subprocess gg qq, is : 

x{K'7("fc - 9, q)T>',tpt{k)V,l?^f,{k + q, -fc, -q)T>'i]^f; {q) -f (A - B)} (A2.3) 

The most difficult problem here is to calculate the trace over spinor and vector indices. The complete expression is 
very unwieldy, and we will therefore satisfy ourselves with 

Trace = 12 x I&p^sxaxbII + 0{pt/ ^/Ts)\, 

where the factor 12 comes from the color algebra. 

As previously, we have three contributions to this process, 

=af (f;,H)+af)(H,H). (A2.4) 

The first term corresponds to the factorization of the parton's structure functions and the "hard" cross-section at a 
given scale, 



da'i^\Q,Q) ialpl dxAdxB 



Pt Pt 



dPtdy 2s3 7o x\x% ^^^^ y/s^ )^{xb ) x 

xGixA,Ql)GixB,Ql)S[xAXB - ^{xAe-y + XBey)]. (A2.5) 
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This term, explicitly dependent on the factorization scale, is kinematically suppressed by three powers of s~^. Two 
of these come from the kinematics of intermediate gluon, and the third from the phase volume which is confined to a 
line in the (z^i, a;B)-plane. 



The second term in Eq. ( A2.4) describes the interaction of a parton with the source, 



2^2 



dptdy 



167r2 



dXA / dxBd{xA 
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ey)0{xB - 



Pt 



hdkt ,GixB,Ql)G'ixA,k^) 



{A ^ B)] 



(A2.6) 



where we have introduced the notation, 



e = {k+ - p+){k- - p-) = s{xA - ^ey){xB " ^e-y). 



This second term still keeps dependence on the factorization scale. Because of finite virtuality of one of the incoming 
fields, the phase volume of th e process is larger, and consequently one of the powers disappears. 
The last term in Eq. (A2.4) accounts for the interaction of two sources, 



dPtdy 
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dx, 



dxBOixA - ^ey)e{xB ~ ^e-y) X 



4:S{lt,Pt,0 



^M'ixA, {It + it)yi)G'{xB, (It - ft)V4) 



{sx\x'^ 



(A2.7) 



Here, both incoming fields are virtual, and the s-channel propagators are smoothed by the total transverse momentum 
of initial fields. The internal integral over d^f represents the two-dimensional Fourier transform of the product of the 
densities (in coordinate space) of two sources. Hence, this term is proportional to the degree of geometrical overlap 
between the colliding nuclei. It does not depend on a factorization scale, and is expected to be the dominant term. 

All three terms do not exhibit any problems at low pt , but a reliable knowledge of the low x behavior of the sources 
may be important for quantitative computations. The presence of a factor p^, which is purely kinematic in its origin, 
guarantees that at high pt the first order differential cross-section will be larger than the Born cross-section in this 
region. 
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FIG. 1. Born diagrams for the one quark and one gluon production. The bold cross labels the line corresponding to the 

"detected" particle with momentum p. The dashed line crosses the Greenians representing densities of the initial (bold) or 
final (thin) states. Numbers near the vertices indicate the type of ordering in the Greenians. The processes: (a) qg q; (b) 
99 9; (a) qq 9- 

FIG. 2. Self-oncrgy-typo first order diagrams for the one quark production. Notation the same as in Fig.l. Arrows label the 
retarded and advanced propagators and show the latest time. 



FIG. 3. Four topologies of the vertex-type diagrams for the one quark production. 



FIG. 4. Twelve vertex-type diagrams of for the one quark production in the first order. 



FIG. 5. Two infrared-divergent "box" diagrams which are not generated by our perturbation theory. 
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